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The dynamical structure factor of the 5* = 1/2 bond-alternating spin chain with a next-nearest- 
neighbor interaction in magnetic field is investigated using the continued fraction method based 
^^ ■ on the Lanczos algorithm. When the plateau exists on the magnetization curve, the longitudinal 

>--^ ' dynamical structure factor shows a large intensity with a periodic dispersion relation, while the 

transverse one shows a large intensity with an almost dispersionless mode. The periodicity and the 
amplitude of the dispersion relation in the longitudinal dynamical structure factor are sensitive to 
the coupling constants. The dynamical structure factor of the S = 1/2 two-leg ladder in magnetic 
field is also calculated in the strong interchain-coupling regime. The dynamical structure factor 
shows gapless or gapful behavior depending on the wave vector along the rung. 
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I. INTRODUCTION 



One-dimensional quantum spin systems have attracted great attention both theoretically and experimentally for 
a long time. Recently, Tonegawa et al. discovered numerically that the S = 1/2 bond-alternating spin chain with 
a next-nearest-neighbor (NNN) interaction has a parameter region where the plateau exists on the magnetization 
curve at half of the saturation value. |jll| The wave function of the lowest-energy state in the plateau region was 
shown to have the periodicity n = 4, ITT indicating that spontaneous translational-symmetry breaking, because the 
periodicity of the Hamiltonian is n = 2. This result is consistent with the necessary condition for the appearance 
'■^ . of the magnetization plateau at half of the saturation value in the 5* = 1/2 spin chain 0, which indicates that the 
^ ' periodicity of the lowest-energy state must be a multiple of 4. On the basis of bosonization technique combined 
O ' with the renormalization-group theory, the magnetization process of this system was investigated in a viewpoint of 
^1 the stability of the magnetization plateau. M A simple explanation for the two- fold degenerate ground-state in the 

■ plateau region was also presented. [g| 
~~^ I The dynamical structure factor (DSF) of quantum spin systems provides us with important information about 

K*" ' the intensity of the magnetic response as a function of energy and momentum transfer. The DSF of 5 = 1/2 one- 
dimensional quantum spin systems in magnetic fields was measured by inelastic neutron-scattering experiments. [0, O 
-^^ ■ The experimental results were compared with numerical results obtained by the exact diagonalization method for the 
r — I finite system, 0, R, 0] and it has been argued that both results are in good agreement. 

f^ The DSF of 5 = 1/2 one-dimensional spin-gapped systems in magnetic fields was investigated by numerical calcu- 

f^ ^ lation P, irn, O. Using the recursion method, the longitudinal DSF of the S = 1/2 Heisenberg chain with a NNN 
^D i interaction was calculated at the half of the saturated magnetization, and the characteristic distribution of the inten- 
sity has been discussed |^. The longitudinal DSF of the S = 1/2 bond-alternating spin chain with a NNN interaction 
was calculated using the numerical diagonalization method llfl, nU. When some adequate parameters are used, this 
model is considered to be an effective model Hamiltonian for CuGeOa. [|l2[ |l^ It has been shown that the large 
intensity lies at the soft mode, pO| , pd| and the distribution of the intensity is remarkably affected by the spin-phonon 
2 ' coupling. [|ll| Using bosonization technique, the DSF of the S = 1/2 one-dimensional spin-gapped systems such as 
Q ' the Heisenberg chain with a NNN interaction, the bond-alternating chain and the two-leg ladder iQ was calculated 
O ■ systematically in magnetic fields [|5|. On the basis of the results, Chitra and Giamarchi have concluded that the 
^ ' features of 5 = 1/2 one-dimensional spin-gapped systems appear not in thermodynamic quantities but in critical and 
• ^H , dynamical properties of the spin correlation function, when the energy gap is collapsed by magnetic fields. p5| 
/\^ '■ In this paper, we calculate the DSF of the S — 1/2 bond-alternating spin chain with a NNN interaction in magnetic 

H , fields, using the continued fraction method based on the Lanczos algorithm. |l6| This model can be regarded as the 
■ - - ' S = 1/2 two-leg spin ladder with a diagonal interaction. In §11, we summarize the numerical method to make this 



paper self-contained. In §[H, we show the DSF of the S = 1/2 bond-alternating spin chain with a NNN interaction, 
turning our attention on the behavior at the magnetization plateau. In § |rv| , the results for the DSF of the S — 1/2 
two-leg spin ladder and the 5* = 1/2 two-leg spin ladder with a diagonal interaction are shown at half of the saturated 
magnetization. We summarize the paper in 



Present address: Research and Development Department, Tochigi Nikon Corporation, 770 Midori, Ohtawara 324-8625. 



II. MODEL AND NUMERICAL METHOD 

Let us consider the following model Hamiltonian for the S — 1/2 bond-alternating spin chain with a NNN interaction 
in magnetic fields, 

H - Ho + H,, (1) 

3 3 1,3 

H, = ~MH, (3) 

where Si.j is the S — 1/2 spin operator on the site j of the sublattice ^Z = 1 or 2), Ji and J2 are the alternating nearest- 
neighbor interactions, and J3 is the next-nearest-neighbor interaction. The chain is divided into two sublattices, which 
consist of every other sites. The distance between neighboring spins is set equal to unity. The magnetization is given 
by M = ^; ■S'f . We set Ji, J2 and J3 > 0. The periodic boundary condition is apphed. Note that the system 
described by the model Hamiltonian (|l|) can be regarded as the S* = 1/2 two-leg spin ladder with a diagonal interaction, 
where Ji, J2 and J3 are the interchain, diagonal and intrachain interactions, respectively. 

To calculate the DSF, we first consider the ladder geometry. In this case, the DSF is defined as, 

n 

= -- lim Im(*o|5'' „ ^ ^ 5*;" „ |*o), (5) 

where z = cj + ie, ^„ is the eigenfunction of Ti, whose eigenvalue is En, vI/q is the ground state eigenfunction with the 
eigenvalue Eq , and qx and qy are the wave vectors along the chain and the rung, respectively. The Fourier transform 
of the spin operator S^^^ is given by 

^t,qy = 7^ XI exP[i(fc j + qyl)]S'^j , (6) 

where fi — x,y, z ot +, — , z. 

We calculate ^0 numerically in the subspace with a given magnetization. The DSF can be expressed 
by the continued fraction in the tridiagonal basis of the Hamiltonian with the initial vector \(f)o) = 



^,^..,J*o)/V(*o|5t.,/5,^.,,J*o);[y 



1 .. {^o\Si:^'S^^,J^o 



S'^{qx,qy,u;)^^- hm ' "' ^^'"'^ '^;'l " , (7) 
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where ai and /3i+i {i = 0, 1,2,- • •) are the diagonal (i -I- l,i + 1) and the sub-diagonal (i -f 2,i + 1) elements of 
the tridiagonal Hamiltonian. The matrix elements ai and (3i are calculated by the Lanczos method. We evaluate 
the continued fraction by the modified Lenz method |lj]. This method enables us to estimate the convergency 
quantitatively. A small but finite e is used instead of taking the limit e ^ 0. We set e = 3.0 x 10^^. 

In magnetic field along the z axis, rotational symmetry around the x and y axes is broken, while rotational 
symmetry around the z axis remains. Hence, S^{qx,qy,Lo) — Sy{qx,qy,oj) ^ S^{qx,qy,uj). We calculate S^{qx,qy,uj) 
and S^{qx,qy,uj) to obtain S'-^{qx:qy,uj) by using S'^{qx,qy,uj) = {l/A)[S'^{qx,qy,Lu) + S^{qx,qy,uj)]. We now define 
the magnetic field where the magnetization changes from Af to M -f 1 or from ill — 1 to M as 

(8) 
(9) 

where El{N, M) is the lowest energy of Tip with N spins and the magnetization M . When we calculate the DSF on 
the magnetization plateau, we have to fix the magnetic field in the plateau region. We use Hi for S~ {qx,qy,i^) and 



:l{n,m + i)-El{n,m) = 


= Hi 


:l(n,m)-El{n,m-i) -- 


= H2 



H2 for S~^ {qx , Qy , to) , indicating that S~{qx,qy,i-o) and S~^{qxTqy,uj) are calculated at the highest and lowest ends of 
the magnetic field on the magnetization plateau, respectively. Then, we fix the magnetic field Hp in H2 < Hp < Hi . 
Shifting S^{qx,qy,uj) and S'^{qx,qy,u!) along the ui axis by Hi — Hp and Hp — i/2, respectively, and adding up 
them thus obtained, we can calculate S^{qx,qy,uj) {= {l/4)[S^{qx,qy,tLi) + S~{qx,qy,uj)]) in a given magnetic field 
Hp. When the magnetization plateau does not appear, the finite size scaling based on conformal field theory is 
worthwhile to calculate the precise value of the magnetic field in the thermodynamic limit for a given magnetization 
m = N/M. pSJ Nevertheless, we take the average of Hi and H2 to determine the magnetic field for a given m, because 
the error in the magnetic field thus obtained shifts slightly the DSF along the lu axis. 

For the one-dimensional geometry, we do not have to distinguish qy from qx. Therefore, (0) and (g) are simply 
described as 

S^{q,L0) = Y.\{^rr\S'i\^^)?5{u ^ [Er.- E^)), (10) 

n 

^'i = ^Ee^P[i(2j + 0'z]^[:r (11) 



On the basis of ( pX]| ) and (|l^) , the calculation of the DSF for the one-dimensional geometry is performed in the same 
way as done for the ladder geometry. 

III. DSF OF THE S = 1/2 BOND-ALTERNATING SPIN CHAIN WITH A NNN INTERACTION 

The DSF of the S = 1/2 bond-alternating spin chains with a NNN interaction is shown in Fig. |l| for to = 1/4 
with {N,M) — (28,7). Note that all figures in this paper have the same scale concerning the intensity of the DSF. 
Parameters are (a) Ji = 1.2, J2 = 0.8, J3 = 0.4, (b) Ji = 1.5, J2 = 0.5, J3 = 0.3, and (c) Ji = 1.1, J2 = 0.9, J3 = 0.3. 
According to ref. 1, the plateau exists on the magnetization curve in Figs. Il|(a) and (b), but do not exist in Fig. 13(c). 
The magnetic fields for m = 1/4 are (a) H = 1.80, (b) H = 1.75, and (c) H ^ 1.751, respectively. In Figs. |l|(a) and 
(b), two divergency of S^{q, lo) in u —> are seen at q = and q — 7r/2, whereas in Fig. ^c) a divergence of S^{q, lo) 
in w — > is seen only at q — 0. The results for S^{q, uj) indicate that in the systems of Figs. |l|(a) and (b), translational 
symmetry is broken and the periodicity is twice as large as that in Fig. yj(c). Thus, the condition for the appearance of 
the magnetization plateau pi is satisfied. The peaks of S^{q,uj) around uj ^ 1.0 in Fig. y(a) and uj ^ 1.3 in Fig. MJd) 
are almost dispersionless. On the contrary, the results for S^{q,uj) change qualitatively depending on the parameters 
as shown in Figs. H(a) and (b). In Fig. 0(a), there are small peaks around uj ^ 1.7. The dispersion relation exhibits 
cosine-like q dependence with the periodicity n and the amplitude ^ 0.7. In Fig. |l|(b), the peaks of S^{q,uj) exhibit 
sine-like q dependence with the periodicity 7r/2 and the small amplitude ~ 0.1. 

We compare the DSF in the absence of the magnetization plateau shown in Fig. |l|(c) with that of the S — 1/2 
Heisenberg chain (Ji = J2 = 1-0, J3 = 0) for m = 1/4 {{N,M) = (28,7)) shown in Fig. |. The magnetic field for 
m = 1/4 in the latter system is evaluated as H ^ 1.584, which is close to the exact value ~ 1.597 [0,^|. The dotted 
lines denote the lower bound of the excitation continuum obtained from the Bcthe ansatz solution |20[] . The results 
for S^{q,u!) and S^{q,uj) in both systems are qualitatively the same, although the additional NNN interaction exists 
in the system shown in Fig. |l|(c). 

IV. DSF OF THE S = 1/2 TWO-LEG SPIN LADDERS 

We show the DSF of the S* == 1/2 two-leg spin ladder system for m = 1/4 with (iV, M) = (28, 7) in Fig. |. Parameters 
are Ji = 2.0, J2 — and J3 = 1.0. Thus, the system is in the strong interchain-coupling regime. The wave vector 
along the rung is fixed to be Ow = and tt. The magnetic field corresponding to rn = 1/4 is evaluated as H ^ 2.768 
from eqs. (H) and (^). In Fig. y[ S^{qx,0, u)) and S^{qx, n, w) are gapless, while S^{qx, vr, uj) and S^{qx,0, uj) are gapful. 
This gapless and gapful behavior is consistent with the asymptotics of the spin correlation functions obtained by the 
hard-core boson representation of the strong-coupling two-leg ladder ||l4| : algebraically decaying ones, (S^qS^q) and 
(Sf^^Sj^^), and exponentially decaying ones, {Sf.^S^.„) and (5'fo5'fo). |24[| 

The strong-coupling S — 1/2 two-leg spin ladder in magnetic field can be mapped onto the S — 1/2 XX Z chain with 
the anisotropy A = 1/2 in effective magnetic field [|l^ ^ ^ |23|. As shown in Fig. ^ (a), the intensity of S^{qx, 0,uj) 
lies almost between the boundaries uj = (3-^/3^3/4)1 sing | and uj = {^V3Js/2)\ sing/2|, which are the lower and the 
upper bounds of the excitation continuum of the S =^ 1/2 XX Z chain with A = 1/2. [ESl In S^{qx,TT,uj), the intensity 
becomes larger along the lower boundary as qx —>■ tt. The magnetic field corresponding to to = 1/4 is given as 
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FIG. 1: The DSF of the S = 1/2 bond-alternating spin chains with a NNN interaction in magnetic fields is shown for m = 1/4. 
Parameters are (a) Ji = 1.2, J2 = 0.8, J3 = 0.4, (b) Ji = 1.5, J2 = 0.5, J3 = 0.3, and (c) Ji = 1.1, J2 = 0.9, J3 = 0.3. The 
magnetic fields for m — 1/4 are (a) H = 1.80, (b) H = 1.75, and (c) H ^ 1.751. In (a) and (b), the plateau exists on the 
magnetization curve. 



Hi/4 = Ji + J3/2 = 2.5 from the effective S == 1/2 XX Z chain with A = 1/2 ||2§. The value ffi/4 = 2.5 is close to 
our numerical estimation H ~ 2.768. 

In S^{qx,T^Tijj) and S^{qx,Q^^^)^ the lower boundaries in the uj — q^ plane show convexly downward and take the 
minima at q^ — 7r/2. The largest intensity of S'^(q3;,7r, w) appears at qx = tt, while that of ^^(ga;, 0, w) appears at 
qx — ^. The origin of the largest intensity can be explained by the fact that the two modes intersect at g^; = tt in 
S^{qx,T^, w), and those intersect at g^; = in S^{qx, 0, lu) as seen in Fig. ^ . These results are qualitatively the same as 
those obtained in ref. 14, where the DSF of the weak-coupling S = 1/2 two-leg spin ladder in < to < 1/2 |2^ was 
calculated analytically using bosonization technique. According to ref. 14, furthermore, it is argued that these largest 
intensities in the gapful DSF come out at (qx,Lu) = {'k^H) in S^{qx,TT,uj) and at {qx^uj) = {0,H) in 5^((72;,0, w). 
From our numerical results, we evaluate that lo ^ H{^ 2.768) at both points in agreement with the results in ref. 
14. It is interesting that the bosonization analysis for the weak-coupling regime well describes the feature of the 
strong-coupling regime. 

Finally, we show the DSF of the S — 1/2 two-leg ladder with a diagonal interaction in comparison with that of the 
S = 1/2 bond-alternating spin chain with a NNN interaction. Used parameters are Ji = 1.2, J2 = 0.8, and J3 = 0.4, 
which are the same as those in Fig. |l|(a). We set m = 1/4 with (TV, M) = (28, 7). The wave vector g^; — tt in the two- 
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FIG. 2: The DSF of the 5 = 1/2 Heisenberg chain in magnetic fields is shown for m — 1/4. Parameters are Ji = J2 ~ 1-0 and 
J3 = 0. 
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FIG. 3: The DSF of the S = 1/2 two-leg spin ladder with (a) qy — and (b) g^ = tt is shown for m = 1/4. The corresponding 
magnetic field is i/ ~ 2.768. Parameters are Ji = 2.0, J2 = 0, and J3 = 1.0. 



leg ladder geometry corresponds to the wave vector q = -k 12 in the one-dimensional geometry. In Fig. ^, two divergent 
peaks appear at q^; = and tt in w — *■ in S'^(ga;, 0,a;), indicating broken translational symmetry. In S^{qx,TT,i^), the 
small dispersionless peaks appear around lu ^ 1.0 and 2.0 as also seen in Fig. |^(a), while in S^ {qx , tt , 1-^) , the peaks 
appear in tj < 1. Since there are no peaks around cj ^ 1.0 and 2.0 in S^{qx, 0, cj), and in w < 1 in S'^{qx, 0, cj), these 
peaks in Fig. 0(b) are governed by the interchain mode qy — ir. This is confirmed by the observation that the intensity 
corresponding to these peaks becomes so weak as q approaches in Fig. 0(a). 
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FIG. 4; The DSF of the S — 1/2 two-leg spin ladder with a diagonal interaction with 
m — 1/4. Parameters are Ji = 1.2, J2 = 0.8, and J3 — 0.4, and the magnetic field is H 



^) <ly = and (b) qy — tt is shown for 
= 1.80. 



V. SUMMARY 

We have calculated the DSF of the S — 1/2 bond-alternating spin chain with a NNN interaction in magnetic fields, 
using the continued fraction method based on the Lanczos algorithm. When the plateau exists on the magnetization 
curve, large intensity with a periodic dispersion relation comes out in ^^(g, oj), whose periodicity and amplitude are 
sensitive to the coupling constants. On the other hand, large intensity with an almost dispersionless mode comes out 
in S^{q, co). When the plateau does not exist, the DSF shows qualitatively the same behavior as that of the S = 1/2 
Heisenberg chain. 

The DSF of the S — 1/2 two-leg ladder for m = 1/4 is also calculated in the strong interchain-coupling regime. 
The DSF shows gapless or gapful behavior depending on the wave vector qy. The gapless modes in S^{qx,0,uj) and 
S'^ {qx , TT , uj) are well described by the S ~ 1/2 XX Z chain with A = 1/2. The gapful behavior in S'^ {qx , tt ^ oj) and 
S'^Cfej 0, a;) is consistent with the results based on bosonization technique for the weak-coupling regime. 
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